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Recent CERN experiments with antiprotonic atoms 
' have renewed interest in the question of the comparison 
^ of the proton and neutron density distribution at the 
surface of large nuclei. It has been known for years that 
hadronic, in particular kaonic atoms, provide a way to 
study the extreme tail of the nuclear density distribution, 
y—( its isospin structure and nuclear correlations there 

■ Two methods have been used , each of which gives 
0^ ' information in the region roughly 2.0 fm beyond the half 
(— I , density radius. 

• 1) An observation of the X-ray cascade in hadronic 
'—I ' atoms and extraction of the atomic level shifts and 
^ widths. These level widths, related to the nuclear ab- 
(— ( sorption of hadrons may determine high moments of the 
• • nuclear density distributions. Information obtained in 
. 1^ this way is rather limited, however. For a given atom 
^ , only one, and in some special cases two level widths and 
^ ' one level shift, can be measured. The latter is in general 
difficult to interpret and provides a check on the hadron- 
nucleon interaction models. 

2) Studies of the nuclear absorption products in par- 
ticular detection of the emitted mesons. In this way one 
can discriminate captures on protons from captures on 
neutrons. In principle, more nuclear information is of- 
fered but it is also more difficult to interpret. The reason 
is that the initial states of capture are not known directly 
and the final mesons may exchange charge or become ab- 
sorbed. 
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The new experiments of Ref. , involving antiprotons, 
have two very strong advantages. Firstly they differ- 
entiate reasonably clearly between the pp annihilations 
and pn annihilations, thus ensuring that the neutron and 
proton distributions can be separately estimated. Sec- 
ondly, they involve antiproton absorbtion very much fur- 
ther from the nuclear half density radius than the earlier 
experiments mentioned above. We shall show that the 
absorption occurs around a region 3.0 fm beyond the half 
density radius. 

The idea of the experiments of Rcf. is to detect 
"cold" nuclei, following the absorption of antiprotons 
in antiprotonic atoms, by radiochemical methods. By 
"cold" here we mean nuclear states of very low energy, 
less than the neutron emission threshold. Characterising 
a nucleus of TV neutrons and Z protons by (TV, Z) the 
reactions are 

p+{N,Z) — > {N,Z -1) + mesons (A) 
p+{N,Z) — > {N -1,Z) + mesons (B) 

Thus (A) involves predominantly interactions of the 
p with protons and (B) similarly with neutrons. It is 
this feature which is the first advantage of this method 
mentioned above. 

The extreme surface nature of the process arises firstly 
from the high orbital angular momentum of the p, but 
this is well known and is exploited in the earlier experi- 
ments j^^. However the new feature, giving increased 
surface emphasis arises from the mesons in (A) and (B). 
On average there are 4-5 mesons emitted and to leave a 
"cold" nucleus they must all avoid hitting the nucleus. 
This can only be achieved in the far surface. 

The basic NN interactions required for this study are 
limited by several phenomenological parameters - range 
of the NN annihilation, absorptive parts of the scattering 
amplitudes, pion production multiplicities, pion momen- 
tum distributions. These are taken from other experi- 
ments while effects of uncertainties must be quantified. 
Once the final and initial states are under control one can 
interpret the (iV —1)/{Z—1) ratios in terms of "neutron 
haloes" or "neutron skins " and attribute quantitative 
meaning to these terms. Qualitatively, the CERN exper- 
iments discover large neutron excess at the nuclear 
surface of several heavy nuclei. It complements similar 
findings in the subcoulomb neutron pickup reactions . 

One purpose of this paper is to provide a description of 
the nuclear capture of the atomic antiprotons. Of main 
interest is that fraction of single nucleon captures which 
leaves cold nuclei in the final states. This is the sce- 
nario of the recent experiment ^ . Difficulties in the way 
are twofold. First, the NN annihilation is a complicated 
process with many degrees of freedom involved. Second, 
the initial atomic state of the antiproton and the final 
states of the residual nuclei are not certain. Fortunately, 
the first difficulty is moderated by a large energy release 
in the p absorption. Hence, a closure over nuclear final 
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states and high energy approximations for the annihila- 
tion mesons may be apphed and these yield classical for- 
mulas for the absorption rates. The latter are expressed 
by integrals of nuclear densities weighted by a probability 
to find the antiproton inside nucleus and a probability to 
find the final nucleus left undestroyed by the annihila- 
tion products. The last two probabilities are calculated 
on the basis of the antiproton and pion optical potentials. 
Their dependence on the initial atomic state, final state 
interactions and parameters of the p and pionic nuclear 
optical potentials are studied in Section II. 

Nuclear models of various complication: Fermi gas, 
shell model, Hartree-Fock and Hartree-Fock-Bogolyubow 
methods are used in Section III to find the nuclear density 
distributions. Comparison with the experiment is done 
and dependence on nucleon binding energies. Coulomb 
barriers, shell structure and correlations is studied. The 
advantages and limitations of the experimental method 
are indicated. 



II. NUCLEAR ABSORPTION OF ATOMIC 
ANTIPROTONS 



Antiprotons bound into atomic orbits cascade down 
to be ultimately absorbed by the nucleus. The latter 
happens at the extreme nuclear surface and the absorp- 
tion probability is significant even at distances as large as 
twice the nuclear radius, Two effects create such 

a situation. First, the free path of antiprotons in nuclear 
matter is less than 1 fm and second, the atomic cascade 
tends to populate states of high centrifugal barrier I. The 
peripherality of capture allows for standard low density 
simplifications: quasi-free scattering and single-particle 
picture of the nucleus. Also, it facilitates the description 
of the final mesons, a vital question for understanding 
the absorption experiments. On the other hand the dis- 
advantage and difficulty inherent in the surface studies is 
related to its sensitivity to range effects. 

This section presents a description of the antiproton 
absorption mechanism. First, a simple phenomenologi- 
cal picture based on the optical potential model is pre- 
sented. Next, two special questions: final state inter- 
actions and range effects are discussed, again in a phe- 
nomenological way. A more detailed justification of the 
phenomenological approach, basic assumptions and lim- 
itations are given in two consecutive subsections. These 
contain rather technical discussion which may be omit- 
ted by readers more interested in the nuclear structure 
results. 

The tool to describe the antiprotonic atomic level shifts 
and widths is an optical potential 1/°p* . The simplest one 
is assumed usually |^,| 1 1 1 in the form 



yopt(j^) = -^t^^piK) (1) 

where /ijviv is the reduced mass, p(R) is a nuclear density 
at a radius R and tj^fj is a complex scattering length. 
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As there is some finite range in the NN interaction the 
density p(R) involved in Eq.(jl|) is not the "bare" nucleon 
density poCR) but a folded one 

p(R) = J dupo(R-u)w(u) (2) 

where u is a formfactor that represents the NN force 
range. For the absorptive part of y°P* the annihilation 
range of 1 fm might be expected from models of the NN 
annihilation but the range in the real part is more difficult 
to control. 

The length ijyTV Eq.(^ is extracted from antipro- 
tonic atomic data. The most precise X-ray measurements 
have been done for the 3d and 4f states in the oxygen iso- 
top es, flpi , and fits to these give tj^j^ of about — 1.5— i2.5 



fm, 1 11 T^. This value yields a deep and strongly absorp- 
tive potential well. At the nuclear centre lmV°^^ would 
be 200 MeV strong and the related free path length would 
be well below 1 fm. However, it should be kept in mind 
that both the form and the strength of y°P' is tested only 
in the surface region. In particular, Imy°P* is determined 
by the atomic level widths, via 

r - A^lmt^j^ /dRp(R) I *^(R) |2 (3) 

where ^jy(R) is the atomic wave function. Since ^'jy ~ 
and only high angular momenta I are available, the 
absorption strength is peaked at the surface. 

The nuclear absorption scenario in ^*Ni is visualised in 
Fig. where some absorption densities W — p \ \^ 
are plotted. There are two special atomic states singled 
out in the capture process. One is the so called "upper" 
level which usually is the last that can be detected be- 
fore the cascading down p is absorbed. One can learn the 
width of this upper state measuring the intensity loss of 
the X-ray transitions. In ^^Ni, and in many other nuclei, 
the nuclear absorption is most likely to happen from this 
level. The next circular state below is called the "lower" 
one. Sometimes, one can measure the shape of the X-ray 
lines feeding this lower state. Such measurements are 
possible when the lower state width is in the range of 
a few keV and additionally the rate of radiation from 
the upper state competes successfully with the upper 
state absorption rate. Chances for antiprotons to reach 
the "upper" levels have not been measured. Presum- 
ably a significant fraction of the antiprotons is absorbed 
in states of much higher principal quantum number n. 
Such a conclusion is reached in the kaonic atoms where 
the X-ray transition intensities for the lower atomic per 
stopped kaon are known to be about 0.2 to 0.6, For 
the studies of n/p ratios however the value of n is not 
very important, one needs to know mainly the distribu- 
tion of capture states in terms of the angular momentum 
I, as at nuclear distances all atomic wave functions of the 
same / differ only by irrelevant normalisations. Calcula- 
tions, |]37|| , indicate the dominance of I = ^uppcr states in 
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the nuclear capture process. We discuss the uncertain- 
ties in this respect in the last section, here we calculate 
the absorption probabilities assuming full occupation of 
the I — ^uppcr + 1 circular level. The capture probabilities 
from several circular states are given in Table I. The lower 
widths are usually larger than upper widths by two or- 
ders of magnitude. That is due to the smaller orbit radii 
and reduced centrifugal barriers. However, the absorp- 
tion density profile is not changed dramatically as may 
be seen from the ratio of these densities given in Fig. 

The ^^Ni nucleus is our reference case. Amongst the 
nuclei tested by the recent p capture experiment, see 
Table II, it is the simplest to describe. 

The localisation of nuclear capture depends on the 
range of NN forces. One way to find this range is to 
use NN potential models, another perhaps more advan- 
tageous, is to fit the atomic and low energy scattering 
data. An early choice was to use charge density profiles 
for the p , [p|j9|JTl[|. This is equivalent to folding a form- 
factor v(u) in Eq.(^) of 0.8 fm rms range. More recently, 
longer ranged gaussian profile formfactors exp(— (r/rp)^) 
have been used, ||l2|,|l3|. Typical best fit values are: rg^ ~ 

1 fm (for ImV) and ror ~ 1-5 fm (for KeV). On the other 
hand, calculations based on the NN potentials yield av- 
erage ranges roi of 0.75 fm up to 1.45 fm, ||lj,|2^, the 
difference beeing due to different handling of the off-shell 
extensions. An effect of the range is shown in Fig. ^. 
A longer NN absorption radius broadens the region of 
nuclear absorption. The related effect on the n/p ratio 
measurement is shown in Fig. |^ and discussed later. 

Optical model calculations based on the TV TV interac- 
tion potentials [p"8|-pTl indicate that the lengths i^jy bear 
no simple relation to the NN S-wave scattering lengths 
which are smaller and repulsive i.e. with real parts pos- 
itive. Thus, Ketj^]^ is of a complicated and uncertain 
structure. At the extreme nuclear surface it reflects a 
long attractive tail of the pion exchange forces, about the 
nuclear radius it may turn to repulsion due to repulsive 
scattering lengths, and is rather uncertain at the nuclear 
matter densities. On the other hand, the phenomenolog- 
ical best fit Imtjy^ represents cumulative effect of the S 
and P wave absorptive amplitudes, and can be well un- 
derstood in terms of the free Imtj^f;. The calculated op- 
tical potentials indicate structure more complicated than 
that given by formula (|^) , but cannot reproduce the data 
as accurately as the latter with the best fit parameters. In 
this calculation we use the phenomenological approach. 

The level widths discussed so far reflect all modes of 
the nuclear absorption of antiprotons. The initial stage, 
an elementary NN annihilation, generates an energy of 

2 GeV of which 3/4 is the kinetic energy taken by the 
final state mesons. The mesons may excite the resid- 
ual nucleus via inelastic scattering and absorption. To 
calculate the total widths one sums over the unobserved 
nuclear excited states. The large energy release and pe- 
ripherality allows to use closure over the nuclear states. 
As a consequence the effective Imt^^ is close to the ab- 
sorptive part of free NN scattering amplitude. That is 
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no longer true when the final nuclear states are limited 
to some particular states, as is the case of radiochemi- 
cal methods that detect "cold" nuclei. The experiments 
in question, allow only final nuclei excited up to the 
neutron separation threshold. Residual nuclei of higher 
excitations would decay by a neutron emission and would 
not be detected by the radiochemistry. In the next sub- 
sections the spectrum of allowed excited states is related 
to the rearrangement - and to interactions of final state 
mesons. 

Now, to explain our aim we give a simplified result, 
which will be proven and refined later. Let s denote 
quantum numbers of the annihilated nucleon (neutron 
or proton), the antiproton atomic orbital (n, l,j) and the 
final state of the residual nucleus (any nucleus or cold nu- 
cleus). We are going to prove a simple expression for the 
partial absorption width Tg corresponding to the "cold" 
nucleus formation. In the zero force range limit this reads 

r, = i^lmt%^ [ d-R I vI/^(R) p p.(R)n(R) (4) 

In this expresion a function Pg is introduced to describe 
formation of the required final states. It is a product of 
two terms Pg = Pmiss * -Pdh- The dominant factor in Pg is 
^miss - the probability that the mesons born at point R 
do not excite the residual nucleus, (missing probability of 
Rcf. [|l|). The other final state factor - Pdh is related to 
the final nucleus rearrangement and happens to be less 
significant. Examples of these functions are given in Fig. 
|. The inte grand in Eq.(^), including the factor, is 
the absorption density for those processes that lead to 
the cold {A — 1) residual nuclei. It is shown as A in Fig. 
^, in comparison to the full absorption density W given 
in Fig. |l| it is shifted to the periphery by almost 1 fm. 

Now, we derive Eq.(^, calculate Pmiss , Ps and study 
the range effects. 

A. Nuclear NN annihilation and final state 
interactions 

The aim of this section is to calculate the rate of nu- 
clear N annihilations that lead to cold final nuclei. This 
is done in several steps: 

1) An amplitude for the NN annihilation into mesons 
tj^ff_,M is assumed and introduced into the nuclear tran- 
sition amplitude in the impulse approximation. 

2) The emission probabilities are calculated and 
summed over the mesonic and nuclear final states. For an 
isolated NN annihilation this procedure would produce 
the absorptive cross section and, via the unitarity condi- 
tion, the absorptive amplitude Imijviv- For nuclear cap- 
tures leading to cold nuclei we limit the summation over 
final states to the states of elastic meson-nucleus scatter- 
ing. This limited summation generates the Imtjviv again, 
but now it is folded over nuclear final state interaction 
factors. 
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3) To simplify our considerations, effects related to 
finite range of the reaction: propagation of the final 
mesonic resonances, recoil effects, nonlocalities due to ex- 
ternal fields and the size of mesonic source are discussed 
at the end of this section. 

Assume, that an antiproton in an n-th atomic state 
annihilates on a nucleon in a single particle state a into 
k mesons with momenta pi , i = I to k. In the impulse 
approximation, the transition amplitude for this process, 
is 



where 5'Ar(x) is the atomic, </5Ar(y) the nuclear, and (fiM 
the mesonic wave function. The latter describes scat- 
tering states and corresponds to the ingoing boundary 
condition. The final state of the nucleus is not specified, 
in the spirit of the impulse approximation it is the ini- 
tial nucleus that is left with a hole in the single particle 
state a. Additional nuclear excitations follow nonelastic 
interactions of the mesons. States generated in this way 
are denoted by index /3 in the mesonic wave functions. 

The transition amplitude tjviv— ►m(^i V-iO known 
in detail. What one needs for atomic studies is, basically, 
the elastic Imt^^ extended off-shell. The momentum ex- 
tension is related to the range dependence expressed in 
terms of the NN relative coordinate x-y. As already 
discussed this is fairly well known, contrary to the range 
dependence in the mesonic coordinates ^. The kinematic 
conditions are special since both the N and N are bound 
and the pair energies fall below the NN threshold. Even 
at these energies, the nuclear momenta reach 1—2 fm~^ 
and the scattering matrix should include at least S and P 
waves with all possible spin states. The relevant partial 
cross sections or partial absorptive amplitudes Imt^jv 
could be calculated from potential models of the NN 
scattering. Such a procedure is adopted in some optical 
potential calculations, even-though the partial 

wave analysis of the NN scattering does not exist. Un- 
fortunately, the problem discussed here is more involved 
and uncertainties are larger. We aim, rather, at a semi- 
phenomenological "effective Imtpffj" as used in the phe- 
nomenological optical potential. 

To calculate the absorption widths, the amplitudes (|^) 
are mod-squared, summed over the final pionic channels 
and integrated over the phase space. One has to sum also 
over the final nuclear states. In this way one arrives at 
an expression for the partial absorption widths Tg 



(5) 




(6) 



where 
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P k '' '' '' 

(7) 

Here, the integration dL means pionic Lorentz invariant 
phase space restricted by the energy conservation and fc 
is the pion multiphcity. For an isolated NN system at 
or below the threshold, P in Eq.(|^) is related by the 
unitarity to absorptive part of the elastic NN scattering 
amplitude 

r - Imi^^(x - y, x' - y')S{R - R') (8) 

where R = (x + y)/2 is the NN CM coordinate. 

For an annihilation inside a nucleus this free space uni- 
tarity relation is no longer true. In the external nuclear 
field the propagation of intermediate particles in Eq.(|^) 
changes. If one is interested in the total annihilation 
rates the summation extends over the mesonic as well 
as over all the nuclear states. In such a case the clo- 
sure approximation applies to nuclear states. It is jus- 
tified by the large energy release, peripherality of cap- 
ture and short annihilation range. Detailed calculations, 
[p"8|-p2|, have been done for N optical potentials with 
models that simulate mesonic channels by complex NN 
potentials. These indicate difficulties of this question at 
nuclear matter densities. The effects of the external nu- 
clear field become noticeable already at p = 0.1/9(0) but 
refiect mainly on the ReF°P*. The antiprotonic atom 
physics discussed here is located safely below this den- 
sity limit. This qualitative picture justifies the success of 
optical potential calculations that can relate the overall 
antiproton absorption rates to the free values of Imi^jy- 
On the other hand, to describe the experimental results 
corresponding to the measurements of Rcf. one limits 
the sum in Eq.(|^) to cold final nuclei. First we discuss 
the chance that the mesons created in the annihilation 
would leave the residual nucleus in such states of low ex- 
citations. Another excitation mode, the rearangement is 
discussed afterwards. 

The spectrum of mesons consists essentially of pious 
correlated in a sizable fraction into p and to resonances. 
These heavy mesons are very broad and after some 1 fm 
propagation range turn to the pions. Their multiplicities 
range from 2 to 8 with an average 4—5 and the average 
momentum is as large as 2 fm~^. Nuclear interactions of 
these pions may be absorptive, inelastic or elastic. Those 
involving a pion absorption occur on two or more nucle- 
ons and produce {A — 2) or lighter nuclei. The inelastic 
processes end up with excited nuclei. These may be the 
{A — 1) nuclei of interest as the dominant mechanisms 
involve single nucleon excitation modes: A and higher 
resonances. For medium and heavy nuclei, and pionic 
energies around A the inelastic cross sections reach 0.5 b 
pO| . The main strengths are located much higher than 
neutron emission thresholds, however. Cross sections for 
excitations of states below these limits are typicaly 1-10 
mb, |pl|. On the other hand the elastic cross sections are 
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very large and reach 1 b, |^^. Hence, the rate for pro- 
duction of cold nuclei is given essentially by the elasticaly 
scattered waves. This allows an optical potential descrip- 
tion. In addition, in the bulk of phase space the pions 
are fast enough to allow an cikonal description. Follow- 
ing this the wave function for each pion is taken in the 
form 

^i7^(pO =exp(ipC-z^(p,e)) (9) 

with S calculated in terms of the pion-nuclcus optical 
potential 

^(P, = dsi^ip^ - C/°pHC + ^s) - p) (10) 

The function S is calculated in a quasi-classical way by 
integrating the local momentum over the stright line tra- 
jectory. Due to nuclear excitations and pion absorptions 
this wave is damped with a rate described by ImS*. The 
latter is generated by absorptive part of the pionic op- 
tical potential Im?7°P*. This damping follows the whole 
path of a pion but the main effect comes from regions 
of large nuclear densities and not the region around the 
birth place ^. We assume that all functions ^(p, £,) are 
related to the central point of annihilation R which is 
the NN CM coordinate. Effects of the source size are 
discussed later, jointly with consequences of the heavy 
mesons propagation range. The assumption on the emis- 
sion of mesons from the central point substantially sim- 
plifies our calculations. With the mesonic wave functions 
which enter Eq.(^, the total momentum of mesons 
P separates to a plane wave form. Some additional de- 
pendence on P is still there but as we show later it is 
rather weak. One consequence is that the NN CM " con- 
servation" 5(R — R') which arises in the free unitarity 
relation (||) is also a good approximation in the nuclear 
case (|^). Now the final state pion interaction factors that 
enter Eq.(0) may be collected into a function 

P'=(R)=<niexp(-5(pi,R))p> (11) 

which is a product of the eikonal factors within each 
multiplicity k. It has to be averaged over the multi- 
plicities and the pionic phase space weighted by some 
unknown momentum dependence generated by t^^v^M- 
The expectation is that the momentum dependence of 
P'^ is weak as compared to the momentum dependence 
of ^ATTv^A/ since the former is determined by the nuclear 
size and the latter by the size of annihilation region. Thus 
one may expected the unitarity condition to hold ap- 
proximately provided the averages of P'' in Eq.(0) are 
factored out. This averaging is now performed and in 
this way one arrives at 

r « Imt^jv(x - y, x' - y')^(R - R')^^miss(R) (12) 
where the "missing probability" is given by 
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P,ni.s(R) = J^Wfe / dLf{p,)Y[ I exp(-5(pi,R)) iV^I^fc / ^^/(?'^) 

(13) 

The integration extends over the restricted Lorentz phase 
space weighted by an experimental multiphcity distribu- 
tion Wk , for k from 2 to 8, |p^ , ^ . In order to check the 
assumptions, some factors f{pi) have been introduced 
into Eg. (|l3|) , while pure phase space and constant t ma- 
trices correspond to / = 1. This probability density 
selects those pionic interactions that do not excite the 
residual nucleus. 

Calculations of Pmiss are performed in a Monte Carlo 
procedure. The optical potential for pions must cover 
a wide momentum range from the threshold up to 0.9 
GeV but the phase space favours a region just above the 
A resonance. This potential is related to the pion nu- 
cleon forward scattering amplitudes and in this way to 
the pion nuclcon cross sections. That method is well es- 
tablished around the A , [|7|. Here, this procedure is 
extended to cover also higher N*i , N^^ resonances which 
are described by Breit-Wigner amplitudes. The two nu- 
cleon absorption mode is taken in a phenomenological 
form Performing these calculations one finds that: 
high energy expansion of the square root in Eq.(|l^) is 
satisfactory, higher resonances cannot be neglected and 
the black sphere limit is a good approximation in dense 
regions. In particular Piniss(R) is changed by less than 10 
percent with an inclusion of the two nucleon absorption. 
The latter beeing of the profile operates in the region 
where the black sphere limit is well fulfiled. An example 
of calculated Piniss(R), given in Fig. ^, is close to a pure 
geometrical estimate that relates it to the solid angle of 
the nucleus viewed from the point R, jl^. Nevertheless, 
the gray zone at the nuclear surface make the effective 
radius of an equivalent black sphere difficult to predict 
off-hand. 

The P,„iss(R) calculated in this way is at best semi- 
quantitative, but the proximity of strong absorption limit 
makes the result fairly independent on the details of the 
annihilation. One question is that the phase space alone 
does not reproduce the experimental momentum distri- 
bution of a single pion, . To remove the discrepancy 
some factors f{pi), which generate the correct distribu- 
tion, have been introduced into average (p^). New, cor- 
rected in this way Pniiss(R) is shown in Fig. ^, it is seen 
to differ only slightly from the pure phase space result. 
As the corrective procedure is uncertain and the change 
is below experimental uncertainties we follow the pure 
phase space averaging. 

In the surface region of interest, the missing probabil- 
ities Pmiss rise linearly with the distance R. That is a 
fortunate result, it makes Pmiss rather insensitive to the 
size and structure of the annihilation region located in 
a small sphere around R. The same applies to effects of 
the p and lo mesons. These may propagate some distance 
to a point R' and decay into pions there. Those events 
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are approximately confined to within a sphere centered 
at R of a radius = velocity * lifetime « 1 fm. Again, the 
linearity of Pi„iss(R.') makes an averaged pionic missing 
probability equal to the Pmiss(R-)- 

Let us turn now to other corrections. The annihila- 
tion happens at the nuclear surface and is confined to a 
region of a small diameter. Nevertheless, effects of the 
nucleus should be considered. These are: the NN cen- 
tre of mass motion, external field, and Pauli principle in 
the intermediate states. At distant surface most of these 
have been found small in the optical model calculations. 
Now we discuss briefly the implications in the A—1 reac- 
tions. The first, elementary effect due to the presence of 
nucleus is the NN centre of mass motion with respect to 
the residual nucleus. It is given by the recoil energy i.e. 
by a Fourier transform of t^^(i? — P'^ /AM) over the the 
NN CM momentum P. It is known in the kaonic atoms 
that a narrow (30 MeV or less) resonance close to 
the threshold would induce a propagation range as large 
as 1 fm and affect strongly the peripherality of nuclear 
capture. Extensive experimental efforts gave no clear ev- 
idence for narrow resonances in the NN system close to 
the threshold. Energy independent transition matrices 
tj^fq^M are assumed here. Hence the annihilation is fast 
and the NN CM may be fixed in the absorption process, 
thus a factor (5(R — R') is justified in Eq.(|l^. On the 
other hand, an effect of the CM motion arises due to de- 
pendence of the pionic functions 5(pi,R) on the total 
momentum of mesons P. It induces some nonlocality in 
the R — R' but the effect enters in a second order of a 
small quantity P/(pion momentum). Numerical studies 
described earlier indicate a nonlocality of a 0.1 fm range, 
negligible in comparison to the 1 fm range effects involved 
in the relative coordinate x — y. 

B. Range effects 

In the limit of zero range NN interactions formula 
for the partial absorption width may be expressed 
in terms of nuclear densities. For finite-ranged interac- 
tions the single particle wave functions involved may be 
reduced only to mixed densities. However, for simplicity 
and historic reasons one wants to have an approximate 
expression in terms of the true densities. At the nuclear 
surface this can be done with a good precision, at least 
for the absorption rates summed over all nucleon states. 
The standard relation, ||2^, that allows it is: 

E ^?(y'X(y) « P.(Y)jo(fcF(Y) I y - y' I) (W) 

a 

where Y = (y + y')/2 , fcp is an effective local Fermi 
momentum which may be calculated in a shell model 
and i — p,n. In a similar way we express the angular 
averaged atomic wave functions ^'5' by 
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2/ + lX!*iV^-^-7VV-;~l -N 



^^(x')nW «l I' {i?"(X - X') + O0L/nB,X'/nB)} 



(15) 



X + Y, 



It is an expansion in the inverse Bohr radius of the orbit 
1/nB and calculation of function W is strightforward. 

Finally, to handle the NN ranges we assume a sepa- 
rable approximation for the scattering matrix in Eq. ( [l^ ) 
tNN = i^(x-y, ro\/2)i^^t;(x'-y',roV2) and gaussian 
form-factors v with some range parameters tq . This al- 
lows for a simple transformation to relative coordinates 
and reduces Eq. (^) into a folded density expression, . 
For the capture rates one has now: 

r. = 4^Imt°^^ / I ^^(Y) |2 viY - X,ro)fxP^mP^,ssi 

H-NN J 

(16) 

where fx is a factor that collects together a large piece 
of finite range effects 

fx= y"duu(u,2ro)jo(fcF(X)u)D(u)du (17) 

With a normalised form-factor v the limit of zero range 
force is fx — 1- For a typical NN absorption range 
rg = 1 fm values of fx ~ 0.5—0.6 are obtained. These are 
almost constant in the nuclear surface region of interest 
and may be taken out of the integral. As our analysis 
involves only ratios of the widths the actual values of fx 
are unimportant. 

The summation in Eq.(|lJ) selects capture events that 
lead to single hole nuclear states. These do not corre- 
spond to the experimental conditions that require cold 
final nuclei or, more precisely, nuclei either in the ground 
states or in states excited below the neutron emission 
threshold r„ . To account for it we limit the sum over 
the initial nucleon states in Eq.([lO|) to those that leave 
the final nucleus with excitation energies less than T„ -I- 2 
MeV where the 2 MeV is allowed for neutron kinetic en- 
ergies. At the surface almost all single particle states of 
sizable overlap with the atomic antiprotons contribute to 
the sum. The experimental cut-off that eliminates deeply 
bound nucleons becomes a (small) correction, which we 
call a "deep hole" factor. Defined as the ratio of lim- 
ited sum of single nucleon densities to the total sum 

Pdh(X) = X^a'* Vq/Sq '/'q and implemented into the 
partial width formula produces 

r^{cold)=4^lmt%j^fx I I vI'jv(Y) |2 r;(Y-X,ro)p.(X)Pdh(X)P„,i3s(^^i^) 

(18) 

This is our result and now we turn to practical, model 
calculations of the basic ingredients in this equation. 
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III. CALCULATIONS OF NUCLEAR DENSITIES 



For the first test we use an asymptotic density (AD) 
model. It follows, essentially, the Bethe-Siemens ap- 
proach p5| ] but it also incorporates larger phenomeno- 
logical input i.e. charge density distribution, neutron 
and proton separation energies and difference of the rms 
radii of proton and neutron densities. At central den- 
sities a Fermi gas of independent protons and neutrons 
is assumed. The Fermi momenta are determined by the 
densities and the Fermi energies are fixed by the separa- 
tion energies. This gives the depth of the potential well 
which in the surface region is extrapolated by the Woods- 
Saxon form. The densities are given by the exponential 
damping of the nucleon wave functions due to the po- 
tential barriers. For protons a Coulomb barrier is added 
and potential parameters (half-density radius c and the 
surface thickness t) are fitted to reproduce the experi- 
mental charge density down to 5 percent of the central 
density. For neutrons the same t is used but c is chosen 
to obtain the rms radius equal (or larger by 0.05-0.10 fm 
in the heaviest nuclei, [^ ) to the proton density rms ra- 
dius. This model is expected to generate average level 
densities, it misses shell effects and correlations. 

As a second method to determine neutron and proton 
densities we have used a selfconsistent Hartree-Fock the- 
ory with the effective two-body Skyrme-type interaction. 
Since our aim in using HF method was rather unusual, 
i.e. to find nucleon densities at the extreme tails of the 
nuclear matter distribution (at distances of 8— 15 fm from 
the center) , a few remarks about its practical implemen- 
tation seem to be in order here. 

The necessary practical condition is the use of a HF 
code not restricting in any way the asymptotic form of 
s.p. wave functions. This condition excludes e.g. all 
codes using the harmonic oscillator basis. In the present 
work we have applied the code solving HF equations on 
the spatial mesh, in which all fields and densities are 
expressed in the coordinate representation. 

The most severe restriction of the presented results is 
the imposed spherical symmetry. It allows enormous sim- 
plification of solution, in particular the HF equation takes 
the form of a differential equation in the radial variable 
for each pair of the conserved s.p. quantum numbers I 
and j. We used 100 mesh points in the radial coordinate, 
in a box of the size of 25 fm and put as a boundary con- 
dition the values of the wave functions at the far end of 
the box equal to zero. 

The density matrix is obtained by summing contribu- 
tions from the lowest s.p. orbits. If necessary, the con- 
tribution of the last orbit is calculated in the filling ap- 
proximation; i.e. an appropriate occupation probability, 
smaller than one, is associated with this orbit. 

The asymptotic form of the radial s.p. wave function 
for large r corresponding to the s.p. energy eigenvalue 

^nl] is 
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Ti-ni]{r) = w;(i)exp(-Kr)/r, (19) 



where n = \J'^\enij\ and wi is a polynomial with 
the dominant term being a constant, so that the con- 
tribution to the asymptotic density is proportional to 
exp(— 2Kr)/r^. For proton orbitals there is an additional 
exponential factor, coming from the Coulomb barrier and 
modifying Kp, which is very important for distances in 
question, i.e., between 7 and 15 fm. Although at very 
large distances from the center the neutron to proton 
density ratio is proportional to exp(2(Kp — Hn)T), where 
K„ and Kp are directly related to neutron and proton 
separation energies, respectively, at distances near the p 
absorption peak usually a few neutron and proton or- 
bitals contribute significantly to the density and a more 
detailed analysis is necessary to evaluate the latter. 

S.p. binding energies (Fermi levels) important in de- 
termination of nuclear density tails are not reproduced 
exactly with existing effective Skyrme forces. In addi- 
tion, calculated spherically symmetric densities for de- 
formed nuclei lack quadrupole correlations which a pri- 
ori may for their own sake distort positions of the Fermi 
levels. Therefore, calculated densities, especially for de- 
formed nuclei, must be treated as approximate, and pos- 
sible sources of error must be kept in mind. In particular, 
the true densities of deformed nuclei may have longer tails 
since their elongated form has to be smeared over Euler 
angles in order to obtain spherically symmetric density 
in the LAB frame. 

In spite of the approximate character of our nuclear 
density calculations we used also HFB theory [ ^J34| in 
order to assess modifications introduced to HF results 
by the residual pairing interaction. In the present case it 
must be distinguished from the very often used HF-f ECS 
method, in which the pairing correlations are included 
using the BCS prescription to self-consistent orbitals. In 
the latter case, the partial occupation of orbitals above 
the Fermi level leads to a non-zero, though usually mi- 
nuscule, occupation of orbitals of positive energy (in con- 
tinuum). Since such orbitals arc not localized this im- 
plies that they dominate nuclar density at large distances 
which is a completely unphysical effect. 

The applied HFB code also uses the coordinate repre- 
sentation and the HFB equations are solved on a spatial 
mesh. The proper analysis of the asymptotic properties 
of two-component quasiparticle wave functions shows js^ j 
that the HFB ground state wave function, even contain- 
ing pairing correlations, is always localized if bound. The 
asymptotic form of the occupied negative-energy quasi- 
particle states is as in the Eq. ( [l9| ) , with k defined in terms 
of the sum Enij—X, where Enij is the quasiparticle energy 
and A is the Fermi energy. 

As the effective force we use the ten-parameter Skyrme 
SkP interaction described in Ref. 1 3^ . It has a virtue that 
the pairing matrix elements are determined by the force 
itself, contrary to other Skyrme-type interactions which 
define only the particle-hole channel. In the paired HFB 
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ground-state the pairing gap is state dependent. As a 
simple pairing gap parameter one can use the pairing 
potential average over the occupied states. 

The "deep hole" corrective factor in the HFB method 
is calculated using the additional condition in the form 
e„ < A, where e„ are the expectation values of the self- 
consistent mean-field matrix (not q.p. energies) and A is 
the HFB Fermi energy. 

The last method we used for determination of nuclear 
densities was chosen to assure correct separation ener- 
gies. A single particle (SP) spherical well, including the 
central and spin-orbit potentials for neutrons and protons 
was assumed. Proper order of s.p. levels is guaranteed 
by the form of the potential. Potential parameters were 
adjusted slightly to obtain the experimental separation 
energies, charge rms radii and, if known, the neutron 
rms radii. 



IV. RESULTS 

In this section we discuss the partial antiproton ab- 
sorption widths: r„(A— 1) and Tp{A— 1) for captures on 
a neutron and proton, respectively, which produce cold 
(A—l) final nuclei. Sum of the two is denoted by T{A—1). 
Experiments determine those partial widths relative to 
the total absorption width r*°*. The data collected in 
Table II consist of two such ratios: o-^-i = T{A— l)/r*°* 
and tT„p = Tn{A — l)/Tp{A — 1). The first one cr^_i is 
a test for a description of the antiproton absorption. In 
particular it checks the weakest points: understanding of 
the final state interactions and knowledge of the initial 
atomic states of capture. If quantitative understanding 
is achieved one can claim control over the region of nu- 
clear surface where the neutron halo is measured. The 
halo itself is seen via the anv ratio. 



A. The (TA-i ratios 

A typical antiproton absorption scenario is visualized 
in Fig. and Fig. ^, which contain some ingredients of 
formula (^q) for the capture widths. The results, given in 
Table I, show cfa-i calculated for some circular atomic 
orbits that are most likely to be the states of nuclear 
capture. The shapes of {A — 1) capture densities are 
determined by the angular momentum Z, strong nuclear 
absorption and Pmiss and thus are rather independent on 
the normalisation of the atomic wave functions i.e. on n. 
Thus, these ratios arc typical to all n states. With the 
angular momentum / increased by one unit, in particular 
from the "lower" to the "upper" and higher I states, the 
increases by about 20%. Thus, the experimental 
data exclude sizable fraction of high I captures but seem 
less restrictive on the states with I lower than /lower , 
where the calculated cr^_i stabilises. On the other hand. 
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cascade calculations done in kaonic, hyperonic and an- 
tiprotonic atoms, show the nuclear capture from 

the latter very low / states to be unlikely. Also, these cal- 
culations indicate an accumulation of the capture prob- 
ability on two or at most three values of I. This result is 
consistent with our capture probabilities, given in Table 
I, and calculated under the extreme assumption that the 
I = ^uppor + 1 circular level is fully occupied at some stage 
of the atomic cascade. On the other hand, calculations 
of Ref. 0] allow a broader distribution with a 20 per- 
cent share of the Supper + 1 and higher I states and a 10 
percent share of Ziowor — 1 values. If that is the real sit- 
uation our results for cta-i and cr„p would rise typically 
by a factor of 1.05. This is the likely uncertainty of the 
calculations in Tables 1, 11 due to poor knowledge of the 
capture orbits. One hopes to clarify some of these points 
by experimental measurements of the cascade intensities 
and absolute cascade intensities in the nuclei of interest, 




The aA-i calculated with the AD and other models are 
consistent with most of the experimental data, shown in 
Table II. This gives some confidence in the validity of the 
final and initial state description. However, there are two 
outstanding discrepancies: Te and Yb. 

Special cases of ^^°Te and ^^^Yb 

The first case is understood qualitatively, the second 
presents a point of specific interest. It is known experi- 
mentally ||26| that a strong E2 mixing, i.e. coupling of the 
atomic and nuclear rotations, occurs in the upper level 
of the ^^°Te atom. It stimulates absorption from the up- 
per level as indicated in Table I, and induces an effect of 
alignments of the nuclear and atomic spins in the states 
admixed to the upper level. Thus the orbital antiproton 
stays closer to the elongated part of the nucleus, as com- 
pared to states of equally averaged orientations. Thus, 
the final pions have a better chance to miss the nucleus. 
Calculations yield some 20% enhancement of the total 
absorption widths due to this effect One expects 

similar enhancement of the cr^-i rate, it is also likely that 
the n/p ratio is higher at the poles of this nucleus. The 
strong E2 mixing happens also in ^^^Yb for high n = 14 
atomic orbits, however. It is not clear as yet what are 
the consequences for the atomic cascade process between 
the n = 14 state and nuclear absorption. It is also not 
clear what correlation of the atomic motion and nuclear 
orientation is induced by this effect. Future experimen- 
tal and theoretical studies |0 will help to elucidate this 
point. 

Antiproton absorption in the heaviest elements, ^'^^Th 
and ^^^U, is accompanied by a nuclear fission of the final 
{A — 1) nuclei which, in principle, may affect the aA-i 
rate. However, in such nuclei the radiative rates domi- 
nate the fission rates for excitations less than the neutron 
emission threshold, js^. In the even-odd nuclei of inter- 
est this domination is even stronger. Thus, the fission 
channel is expected to change the (ta-i only a little, and 
this is apparently borne out by the data in Table 11. 
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B. The CT„p ratios 



The partial absorption widths are proportional to ef- 
fective absorptive amplitudes Imt^^ for the pn and pp 
pairs. These are not well known, although some average 
values follow from the optical potential phenomenology. 
The number required for neutron halo studies is a ratio 
Rnp = Imtlpn) /hnl(pp) which may be taken from other 
experiments. One value Rnp = 0.63 has been obtained 
by Bugg, |3| , from measurements of charged pions emit- 
ted in the p absorption in Carbon. Difficulties arise since 
it includes effects of final state mesonic interactions and 
the inherent uncertainties of the charge exchange reac- 
tions. This value of Rnp generates mild disagreement 
with the data of Ref. 0] for all the nuclei and all nuclear 
models used here. The results given in Table I should be 
compared to the experimental data in Table II, similar 
discrepancies are generated by other models. A different 
result Rnp — 0.81(3) follows from the stopped antiproton 
absorption in deuterium, . This value is free from the 
pionic effects, but the deutron kinematic conditions, in 
particular the binding energy, differ from those met at 
the nuclear surface. Another value of Rnp obtained in 
■^He is smaller and energy dependent. At rest, a num- 
ber Rnp — 0.48(10) has been deduced from rather in- 
volved analyses of the final state mesonic interactions in 
the three nucleon systems [ p2[ . 

In this work we fix Rnp from a best fit to our simplest 
nucleus which is ^^Ni. Our nuclear models yield similar 
results in this case and the fitted Rnp is very close to the 
value obtained from the deuteron. We shall use the latter 
in our analysis. 

The results for CT„p are collected in Table II. It is clear 
that the crudest model of asymptotic density strongly 
overestimates the Pn/Pp ratios at large distances. This 
property has been known already from the neutron 
pickup studies, ||l^. By the same effect, the AD model 
produces too large aA-i in the heaviest nuclei Th and U. 
The physics behind it is quite transparent: 

1) It is vital to have correct separation energies but 
these are not the whole story. 

2) The Coulomb barriers enhance anomalously the 
Pn/Pp ratio at large distances. That has to be off-set 
by shell effects (angular momentum barriers) and corre- 
lations. 

3) Proper setting of the neutron skin defined in terms 
of mean squared radii Rms(neutrons) — Rmsiprotons) 
does not determine the " neutron halo" . The latter are 
understood here in terms of Unp i.e. ratios of high mo- 
ments of density distributions. 

What are the moments involved in the haloes measured 
by the radiochemical experiment? For zero range inter- 
actions and Pniiss ~ 1 these are the "Barret moments" 
i.e. 2/-th moments due to centrifugal barriers corrected 
for the atomic wave functions. The Pmiss and Pdh in- 
crease the order of the moments approximately by two 
units. On the other hand, the annihilation range effects. 
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i.e., the folding, introduces moments smaller by two, four 
and more units. The joint effect is best estimated by the 
dominant 21 density moment involved. As we see from 
Table I these are very high moments ranging from 10 in 
Ni to 18 in U. 

The HF method with the SkP force gives roughly cor- 
rect neutron separation energies for (nearly) spherical 
systems ^^Ni, ^^Zr, ^^''Sm, underbinds the last neutron 
in ^^Ru by 1.3 MeV but generally underestimates proton 
separation energies, e.g. by 1.1 MeV in ^''^Sm and by 3 
MeV in ^*Ni. This statement is qualitatively true also 
for the rest of deformed nuclei, with the one exception of 
^^''Sm where the last neutron is underbound by 2.9 MeV. 
As we have checked in a separate calculation the Skyrme 
force SIII does not improve description of Fermi energies 
in the studied nuclei. 

Comparison between the data and HF results shows 
clear disagreement for Yb and ^"^^Sm and less pronounced 
one for Te, Th, and U nuclei. On the basis of compar- 
ison of calculated and experimental separation energies 
one can expect qualitatively corrections to the calculated 
(7A-1 and cr„p ratios (remembering that like errors in pro- 
ton and neutron Fermi energies tend to compensate each 
other for the anp)- In all presented cases except ^'^''Sm 
and ^^^Yb they go in the right direction. 

In order to better understand the asymptotic HF den- 
sities one can look closer at how many orbitals contribute 
to it at large distances and to test by means of the for- 
mula (14) the sensitivity of the Pn/ Pp ratio to shifts in s.p. 
energies. In ^*Ni, at r = 7 fm, two orbitals give about 
35% of the neutron density each, 3 others give about 10% 
each. For protons, there is one orbital contributing 56% 
(the highest one) and three other contributing about 12% 
each. At r = 15 fm, the last occupied orbitals contribute 
89% and 91% to the neutron and proton densities, re- 
spectively. Of course, the heavier the nucleus, the farther 
the asymptotic region and the more orbitals contribute 
to one-body densities in the range 8-15 fm. In ^''^Sm, 
at r = 7 fm, 4 states contribute more than 10% each to 
the neutron density (the highest percentage being 24%), 
four others contribute more than 5% each. Even at 15 
fm, there are still 4 neutron orbitals contributing signifi- 
cantly to Pn (32%, 24%, 20% and 16%), the reason being 
that their s.p. energies differ not more than 4 MeV. In 
^^^U, at 15 fm, one must account for 4 neutron orbitals 
while only two proton orbitals contribute 79% and 15% 
of the Pp. Clearly, it is not only a distance r but also 
the (sub) shell structure which decides how many orbitals 
contribute. 

In order to estimate the effect of s.p. energies on the 
density ratio Pn/ Pp we have used the following procedure. 
Contributions to the density at r = 8 fm were used as a 
data and the propagation to r = 15 fm was performed us- 
ing Eq. (p^ , with account for the Coulomb barrier. The 
resulting densities for ^^Zr are then larger than the exact 
HF densities roughly by the factor 7/4 and 9/7 for neu- 
trons and protons, respectively. (The Pn/ Pp ratio is then 
161 instead of 115.) This error comes from the influence 
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of the polynomial wi in Eq.(p^ but we concentrate on 
the effect of changing s.p. energies on the so calculated 
Pn/ Pp 8tt r — 15 fm. The decrease of the neutron energies 
by 1 MeV changes this ratio to 94, the decrease of the 
proton energies by 1 MeV rises the ratio to 234 and the 
simultaneous decrease of both proton and neutron ener- 
gies by 1 MeV slightly decreases the ratio to 137. This 
gives some feeling as to the sensitivity of the Pn/ Pp to 
the neutron and proton binding energies. 
Special case of ^^Ni 

This is our reference nucleus. The input data is certain 
and all model calculations produce consistent results. We 
use Rms{n) = 3.734 and Rms{p) = 3.710 , which in 
the SP approach produce aA-i — 1-05 and cr„p = 0.88 . 
That is close to the HF and HFB results from Table II. 
The latter two methods are not perfect, HF underbinds 
protons by 0.5 and neutrons by 0.3 MeV. These errors 
are of no significance, however, since the absorption is 
spread over four neutron and four proton s.p. states. 

Special case of ^''^Sm 

This nucleus displays a proton halo. Qualitatively, it 
might be due to the closed neutron shell. The separation 
energy of neutrons is large (10.6 MeV) with respect to 
a small (6.2 MeV) one for the protons. However, this 
proton halo effect is not reproduced in our calculations. 
Boths results for cr„p given in Table II and plots of the 
neutron / proton density ratios given in Fig. indicate a 
neutron excess at the surface. The separation energies in 
our nuclear models are either fitted to the experimental 
(AD or SP) or well reproduced for neutrons. For protons 
the HF model underestimates the experimental value by 
0.9 MeV and that should even enhance the proton tail 
over the real one. This case indicates again that it is not 
only the separation energy that matters for the nuclear 
tail. The experimental result in ^*'*Sm is not understood 
and opens the case for more exotic speculations. 

One obvious effect of inclusion of pair correlations is a 
change in the Pn/Pp ratio in the tail of the density fol- 
lowing from the change in the Fermi energy. At smaller 
distances however, this ratio may change in the opposite 
direction if other levels than the last one contribute to 
the density. It turns out that in all studied cases the 
change in the Pn/Pp due to pairing is small up to 14 
fm. A much more pronounced pairing effect on both ra- 
tios (T„p and cr^_i comes from the Pdh factor which e.g. 
in ^^^Th changes from 0.52 (no pairing) to 0.69 (with 
pairing) at the total absorption peak at r = 8 fm. This 
change more than balances a decrease in Pn/ Pp ratio pro- 
viding for larger values of (ja-i = 0.109 and cr„p = 4.65 
(see Table II). The same Pdh factor is responsible for an 
increase of both ratios in ^'^^U, to 0.100 and 4.31, re- 
spectively, while the smaller Pdh leads to smaller cta-i 
and i7„p in ^^Zr (0.117 and 2.40, respectively). Pairing 
changes also (ja-i for -'^^'^Sm to 0.106 and (T„p for ^^''Sm 
and ^^^Yb to 2.96 and 3.50, respectively. Changes due 
to pairing are nearly none for other nuclei. 
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V. CONCLUSIONS 



The radiochemical method which detects the products 
of miclcar capture of antiprotons is a valuable source of 
information on the relative neutron /proton density dis- 
tribution on the extreme tail of nuclear surface. The 
main features which we want to stress are: 

1) The nuclear regions tested are more peripheral than 
those studied by the X ray measurements in hadronic 
atoms. One measures essentially the 21 moments of the 
density distributions where I is the angular momentum 
of "upper levels". These moments are as high as 18 in 
the heaviest nuclei. 

2) There are special cases of alignment of nuclear and 
atomic angular momenta formed by the E2 mixing which 
display higher n/p ratios and higher rates of cold single 
nucleon captures. These may test the composition of 
the pole regions in deformed nuclei. Further studies are 
recommended. 

3) The uncertainty in the initial atomic state of capture 
is kept under a fair control by the ua-i rates. However, 
additional X-ray experiments would be helpful to clarify 
this question. 

4) Strong neutron haloes are observed in heavy de- 
formed nuclei. These are not determined by the binding 
energies and Coulomb barriers alone. The shell effects 
(angular momentum) are also important. Few (two or 
three) of the highest nucleon orbitals contribute most to 
the cold capture rates. 

5) An interesting case of a proton halo is found in 
^^''Sm. It is not understood in terms of single particle 
models, and may signal strong nuclear correlations in the 
surface region. 

6) Apart from the ^"^^Sm case the nuclear models re- 
produce the qualitative features of the observed haloes. 

We thank J. Jastrz§bski for his constant interest and 
collaboration and P.E. Hodgson for his advice and for 
arranging the collaboration that was kindly supported 
by The Royal Society and KBN. 



[1] J. Jastrzgbski, H. Daniel, T. von Egidy, A. Grabowska, 
Y.S. Kim, W. Kurccwicz, P. Lubinski, G. Ricpc, W. 
Schmid, A. Stolarz, and S. Wycech, Nucl. Phys. A558, 
405c (1993); P. Lubinski, J. Jastrzgbski, A. Grochulska, 
A. Stolarz, A. Trzciriska, W. Kurcewicz, F.J. Haxtmamn, 
W. Schmid, T. von Egidy, J. Skalski, R. Smolariczuk, S. 
Wycech, D. Hilscher, D. Polster, and H. Rossner, Phys. 
Rev. Lett. 73, 3199 (1994). 

[2] D.H. Wilkinson, Phil. Mag. 4, 215 (1959); P.B. Jones, 
Phil. Mag. 3, 33 (1958). 

[3] W.M. Bugg, G.T. Condo, E.L. Hart, H.O. Cohn, and 
R.D. McCuUoch, Phys. Rev. Lett. 31, 475 (1973). 



20 



[4] D.H. Davis, S.P. Lovell, M. Csejthey-Barth, J. Sacton, G. 
Schorochoff and M. O'Reilly, Nucl. Phys. Bl, 434 (1967). 

[5] E.H.S. Burhop, Nucl. Phys. Bl, 438 (1967); Nucl. Phys. 
B44, 445 (1972). 

[6] C.E. Wiegand and G.L. Godfrey, Phys. Rev. A9, 2284 
(1974); R. Soki and C.E. Wiegand, Ann. Rev. Nucl. Sc. 
25, 241 (1975). 

[7] P. Robcrson, T. King, R. Kunselman, J. Miller, R.J. Pow- 
ers, P.D. Barnes, R.A. Eisenstein, R.B. Sutton, W.C. 
Lam, C.R. Cox, M. Eckhouse, J.R. Kane, A.M. Rush- 
ton, W.F. Vulcan and R.E. Welsh, Phys. Rev. C5, 1954 
(1977). 

[8] Th. Kohler, P. Bliim, G. Biiche, A.D. Hancock, H. Koch, 
A. Kreissl, H. Poth, U. Raich, D. Rochmamn, G. Backen- 
stoss, Ch. Findeisen, J. Repond, L.Tauscher, A. Nilsson, 
S. Carius, M. Suffcrt, S. Charalambus, M. Cliardalas, S. 
Dedoussis, H. Daniel, T. von Egidy, F.J. Hartmann, W. 
Kanert, G. Schmidt, J.J. Reidy, M. Nicholas and A. Wolf, 
Phys. Lett. B176, 327 (1986). 

[9] D. Rohmann et al., Z. Phys. A325, 261 (1986). 
[10] H.J. Korner and J.P. Schiffer, Phys. Rev. Lett. 27, 1457 
(1971). 

[11] C.J. Batty, Nucl. Phys. A372, 433 (1981). 
[12] C.J. Batty, Phys. Lett. B189, 393 (1987). 
[13] E. Friedman and J. Licliteristadt, Nucl. Phys. A455, 573 
(1986). 

[14] The CERN proposal. Experiment PS 209. 
[15] J. Cugnon and J. Vandermeulen, Ann. Phys. (Paris) 14, 
49 (1989). 

[16] P. Arrnanteros ct al, CERN/PSCC/80-101,1980. 

[17] H.R. Clover, R.M. DcVries, N.J. DiGiacomo, and Y. 

Yariv, Phys. Rev. C26, 2138 (1982); J. Roy et al, in 

Proceedings of International Seminar on NN, Syracuse, 

1975, edited by T.E. Kaloreropoulos. 
[18] A.M. Green and S. Wycech, Nucl. Phys. A377, 441 

(1982). 

[19] T. Suzuki and H. Narumi, Nucl. Phys. A426, 413 (1984). 
[20] S. Adachi and H.V. von Geramb, Acta Phys. Aust. 

XXXVII, 627 (1985); Nucl. Phys. A470, 461 (1987). 
[21] O. Dumbrajs, H. Heiselberg, A.S. Jensen, A. Miranda, 

G.C. Oades, and J.M. Richard, Nucl. Phys. A457, 491 

(1986) . 

[22] J. Kronenfeld, A. Gal, and J.M. Eisenberg, Nucl. Phys. 

A430, 413 (1984). 
[23] A.M. Green and S. Wycech, Nucl. Phys. A467, 744 

(1987) . 

[24] X. Campi and A. Bouyssy, Phys. Lett. 73B, 263 (1978). 

[25] H. Bethe and Ph. Siemens, Nucl. Phys. B21, 587 (1970). 

[26] S. Wycech, F.J. Hartmann, H. Daniel, W. Kanert, H.S. 
Plendl, T. von Egidy, J.J. Reidy, M. Nicholas, L.A. Red- 
mond, H. Koch, A. Kreissl, H. Poth, and D. Rohmann, 
Nucl. Phys. A561, 607 (1993). 

[27] L.C. Liu and CM. Sliakin, Phys. Lett. B78, 389 (1978). 

[28] J.N. Ginocchio, Phys. Rev. C17, 195 (1978). 

[29] R. Bizzari, P. Guidoni, F. Marcelija, F. Marzano, E. 
CastelM, M. Sessa, Nuovo Cim. 22A, 225 (1974); T.E. 
Kalogeropoulos, Phys. Rev. D22, 2585 (1980). 

[30] D. Ashery, Nucl. Phys. A355, 385 (1980); D. Ashery, I. 
Navon, G. Azuleos, H.K. Walter, H.J. Pfeifer and F.W. 
Schlepiitz, Phys. Rev. C23, 2173 (1981). 



21 



[31] J. P. Shifler, Nucl. Phys. A355, 339 (1980). 

[32] F. Balestra et.al, Nucl. Phys. A491, 572 (1989). 

[33] N.N. Bogolyubov, JETP (Sov. Phys.) 7, 41 (1958); Sov. 

Phys. Usp. 2, 236 (1959); Usp. Fiz. Nauk 67, 549 (1959). 
[34] P. Ring and P. Schuck, The nuclear many-body problem 

(Springer- Vcrlag, New York 1980). 
[35] J. Dobaczewski, H. Flocard, and J. Treiner, Nucl. Phys. 

A422, 103 (1984). 
[36] G.R. Satchler and W.G. Love, Phys. Reports C55, 183 

(1979). 

[37] Y. Eisonborg and D. Kcsslcr, Phys. Rev. 130, 2352 
(1961); A.D. Martin, Nuovo Cim. 27, 1359 (1963); M. 
Leon and R. Seki, Phys. Rev. Lett. 32, 132 (1973). 

[38] J. A. Wheeler, Handbook of Physics, edited by Condon 
and Odishow, (McGraw-Hill Book Comp., 1958), p.l77. 



22 



TABLE I. Atomic results. Column 2 contains the principal q. number n and the angular momentum /. For the 
remaining columns: c.p. is the nuclear capture probability calculated under the assumption that the circular atomic 
state n — nuppcr + 1 is fully occupied, aA-i= r"*~^/r*°*is the branching ratio for the cold capture and (Jnp— Fn/Fp 
is the ratio of captures on neutrons to protons. The AD model and Rnp = 0.63 were used. 



ELEMENT 


n 


1 


c.p. 




(Jjip 


58Ni 


4 


3 


0. 


0.095 


0.69 




5 


4 


0.16 


0.097 


0.69 




6 


5 


0.83 


0.110 


0.70 




7 


6 


0.01 


0.150 


0.71 




8 


7 


0. 


0.220 


0.71 


96 2r 


6 


5 


0.24 


0.106 


4.67 




7 


6 


0.72 


0.128 


5.30 


i30Te 


7 


6 


0.05 


0.096 


1.77 




8 


7 


0.93 


0.122 


2.00 


Sm 


7 


6 


0.01 


0.075 


1.39 




8 


7 


0.75 


0.085 


1.46 


Sm 


7 


6 


0.01 


0.087 


3.65 




8 


7 


0.75 


0.099 


3.98 


176Yb 


8 


7 


0.23 


0.097 


3.34 




9 


8 


0.75 


0.124 


4.07 


232Th 


7 


6 


0. 


0.073 


3.94 




8 


7 


0. 


0.091 


4.64 




9 


8 


0.31 


0.098 


5.00 




10 9 


0.69 


0.127 


6.20 


238u 


9 


8 


0.29 


0.106 


6.55 




10 9 


0.71 


0.138 


8.24 



TABLE IL Comparison of nuclear models. Experimental and calculated results for aA-i and cr„p are given. Cal- 
culations for atomic orbitals weighted as in Table I are done with i?„p = 0.82 . 



Exp. Ill AD HF HFB 

CA-l CTrtp CTA-I <7np O' A-1 O'np O' A-1 O'np 

58Ni 0.098(8) 0.9(1) 0.11 0.90 0.110 0.785 0.110 0.781 

96Zr 0.161(22) 2.6(3) 0.12 4.9 0.125 2.54 0.117 2.40 

96Ru 0.113(17) 0.8(3) 0.10 1.7 0.099 0.944 0.099 0.955 

130Tc 0.184(36) 4.1(1) 0.12 2.6 0.124 3.14 0.123 3.22 

i44Sm 0.117(20) < .4 0.09 1.9 0.094 1.38 0.092 1.36 

i^'^Sm 0.121(20) 2.0(3) 0.10 5.1 0.110 3.34 0.106 2.96 

"^Yb 0.241(40) 8.10(7) 0.12 4.8 0.111 3.23 0.109 3.50 

232Th 0.095(14) 5.4(8) 0.12 7.6 0.087 3.80 0.109 4.65 

238U 0.114(9) 6.0(8) 0.13 10 0.092 4.09 0.100 4.31 
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FIG. 1. The total antiproton absorption densities from the "upper" n — 6 , I = 5 orbit in ^*Ni: Wi for the NN annihilation 
range ro = 1 fm and Wa for the range ro = 0.75 fm. The dot-dashed line is a ratio of two Wi for subsequent circular n = 5 and 
n — 6 atomic states, po is a "bare" neutron density. Normalisations are arbitrary. 



FIG. 2. The CA — 1) "cold" antiproton absorption density on a neutron from n = 6 circular orbit in ^^Ni. Ai given by the 
integrand of Eq.M) for the NN annihilation range ro = 1 fm and ^4^ for the range ro = 0.75 fm. po is a "bare" neutron density. 
Normalisations are arbitrary. Missing probabilities (left scale): Pmiss continuous is due to phase space alone, Pmiss dash-dotted 
is calculated with corrections for the experimental pion momentum distribution. The flat dashed curve is Pdh from the HFB 
model. 



FIG. 3. The ratios of neutron and proton densities calculated with several nuclear models in the Sm nucleus. Dashed 
line-HF, dotted-SP, dot dashed-AD, continuous line-HFB. The {A — 1) "cold" absorption density peaks at 8.6 fm marked with 
a dot while its bulk is located between 7.2 and 10.2 fm. The experimental anp < 0.4 . 
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